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ABSTRACT 

A new type of wind - a conical wind - has been discovered in axisymmetric magneto- 
hydrodynamic simulations of the disk-magnetosphere interaction in cases where the 
magnetic field of the star is bunched into an X-type configuration. Such a configura- 
tion arises if the effective viscosity of the disk is larger than the effective diffusivity or 
if the accretion rate in the disk is enhanced. Conical outflows flow from the inner edge 
of the disk into a narrow shell with half-opening angle of 30 — 45°. The outflow carries 
<~ 0.1 — 0.3 of the disk mass accretion rate and part of the disk's angular momentum. 
The conical winds are driven by the gradient of the magnetic pressure which exists 
above the disk due to the winding of the stellar magnetic field. Exploratory 3D sim- 
ulations show that conical winds are symmetric about rotation axis of the disk even 
if the magnetic dipole is significantly misaligned with the disk's rotation axis. Coni- 
cal winds appear around stars of different periods. However, in the case of a star in 
the "propeller" regime, an additional - much faster component appears: an axial jet, 
where matter is accelerated up to very high velocities at small distances from the star 
by magnetic pressure force above the surface of the star. The simulations are done in 
dimensionless units and are applicable to a variety of the disk-accreting magnetized 
stars: young stars, white dwarfs, neutron stars, and possibly black holes. For the case 
of young stars, conical winds and axial jets may appear in different cases, including 
Class I young stars, classical T Tauri stars, and EXors. In EXors periods of enhanced ac- 
cretion may lead to the formation of conical winds which correspond to the outflows 
observed from these stars. 
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1 INTRODUCTION 

Jets and winds are observed in many objects from young stars 
to systems with white dwarfs, neutron stars, and black holes 
(e.g., Livio 1997). For example, jets have recently been ob- 
served from the accreting neutron star Circinus X-l (Heinz 
et al. 2007). The jet shows a conical structure with a half- 
opening angle about 30°. This conical structure could be ex- 
plained by precession of the jet. However, 10 years of obser- 
vations show that orientation of the jet has not changed (Tu- 
dose et al. 2008). Jets and outflows are also observed from 
symbiotic stars where matter flows from the vicinity of white 
dwarfs during periods of enhanced accretion (e.g., Sokoloski 
& Kenyon 2003). It is often the case that outflows are as- 
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sociated with enhanced accretion (e.g., Cabrit et al. 1990; 
Lovelace, Romanova, Newman 1994). 

A large body of observations has been accumulated on 
outflows form young stars at different stages of their evolu- 
tion from very young stars where powerful jets are observed 
up to classical T Tauri stars (CTTSs) where weaker jets and 
winds are observed (see review by Ray et al. 2007). A sig- 
nificant number of CTTS show signs of outflows in spectral 
lines, in particular in He I (Edwards et al. 2006; Kwan, Ed- 
wards, & Fischer 2007). High-resolution observations show 
that outflows often have an "onion-skin" structure, with high- 
velocity outflows in the axial region, and lower-velocity out- 
flow at larger distance from the axis (Bacciotti et al. 2000). 
High angular resolution spectra of [Fell] A 1.644/im emission 
line taken along the jets of DG Tau, HL Tau and RW Auriga 
reveal two components: a high- velocity well-collimated ex- 
tended component with velocity v ~ 200 — 400 km/s and a 
low-velocity v ~ 100 km/s uncollimated component which 
is close to the star (Pyo et al. 2003, 2006). High-resolution 
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observations of molecular hydrogen in HL Tau have shown 
that at small distances from the star the flow shows a coni- 
cal structure with outflow velocity ~ 50 — 80 km/s (Takami 
et al. 2007). EXors represent an interesting stage of evolution 
of young stars where enhanced accretion leads to period of 
enhanced outflows (e.g., Brittain et al. 2007) 

Different models have been proposed to explain out- 
flows from CTTSs (see review by Ferreira, Dougados, & Cabrit 
2006). The models include those where the outflow originates 
from a radially distributed disk wind (Konigl & Pudritz 2000; 
Ferreira et al. 2006) or from the innermost region of the ac- 
cretion disk (e.g., Lovelace, Berk & Contopoulos 1991; Livio 
1997). Further, there is the X-wind model (Shu et al. 1994; 
2007; Najita & Shu 1994; Cai et al. 2008) where most of 
the matter outflows from the disk-magnetosphere boundary. 
This work focuses on the outflows from the inner part of the 
disk or disk-magnetosphere boundary. We discovered condi- 
tions which favor conical outflows and performed axisymmet- 
ric and exploratory 3D MHD simulations for both slowly and 
rapidly rotating stars including stars in the propeller regime. 

Earlier we investigated outflows from accreting rapidly 
-rotating magnetized stars in the "propeller" regime (e.g., II- 
larionov & Sunyaev 1975; Lovelace, Romanova & Bisnovatyi- 
Kogan 1999; Alpar & Shaham 1985) where the magneto- 
sphere rotates faster than the inner regions of the disk and the 
magnetosphere transfers its angular momentum to the disk 
matter (e.g. Lovelace et al. 1999). Simulations show that most 
of matter flows from the inner parts of the disk to wide-angle 
open cones. At the same time most of the angular momentum 
and rotational energy of the star outflows along open stellar 
magnetic field lines in a collimated, magnetically dominated 
jet of low-density and high-velocity (Romanova et al. 2005; 
Ustyugova et al. 2006). 

In the present study we were able to obtain conical winds 
in cases of slowly rotating stars, where new disk matter ar- 
rives after period of lower accretion. In this situation the 
poloidal field lines are bunched into an X-type configuration 
and persistent conical winds are observed. Compared to X- 
winds proposed earlier by Shu and collaborators (e.g. Shu et 
al. 1994) conical winds are strongly non-stationary and they 
may originate from a star with any rotation rate. Such winds 
may appear and blow during periods of enhanced accretion, 
or may appear in any situation where the inward accretion 
velocity of the disk is higher than outward diffusion of the 
magnetic field. Simulations were done in dimensionless form 
and such winds may appear around magnetized stars of dif- 
ferent sizes from young stars to neutron stars, and also possi- 
bly around magnetized black holes. 

In §2 we describe the model. In §3 we describe the 
physics of conical outflows. In §4 we discuss dependence of 
the results on different parameters. In §5 we show exploratory 
3D simulations of conical outflows. In §6 we discuss appli- 
cation of the model to EXors and CTTSs. In §7 we compare 
conical winds with propeller-driven winds and other types of 
winds. In §8 we present discussion and conclusions. 



2 MODEL 

We have done axisymmetric and exploratory 3D MHD simu- 
lations of the interaction of an accretion disk with the mag- 
netosphere of a rotating star. The corotation radius of the star 
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Table 1. Sample reference values of the dynamical quantities used in 
our simulations. 



r cr = (GM*/fiJ) 1/3 is in most cases larger than the magne- 
tospheric radius r m which is determined by the balance be- 
tween the pressure of the star's magnetic field and the ram 
pressure of the disk matter. 

The model we use is similar to that described in Ro- 
manova et al. (2005) and Ustyugova et al. (2006). Specifi- 
cally, (1) a spherical coordinate system (r, 0, (f>) is used which 
gives high grid resolution near the dipole; (2) the complete 
set of MHD equations is solved to find the eight variables 
(p, v R , vg, vj,, B r , B g , B^, e) (with e the specific internal 
energy); (3) a Godunov-type numerical method is used; (4) 
an effective viscosity and effective magnetic diffusivity are in- 
corporated the code using an a— model controlled by dimen- 
sionless parameters a v and ay. We varied both parameters in 
the range 0.01 — 1. 

2.1. Initial and Boundary Conditions. Initially we sug- 
gest that the disk is far away from the star (beyond the ex- 
ternal boundary of the simulation region), so that the non- 
perturbed dipole field of the star initially occupies the en- 
tire simulation region and the magnetic field of the star is 
a dipole: B = [3(/x ■ r)r — /xr 2 ]/|r| 5 , where fi is the star's 
magnetic moment. The simulation region is filled with a low- 
density corona and is initially non-rotating. The corona is 
isothermal and in thermal equilibrium, so that the density and 
pressure distributions are: p — p exp[GM/{c 2 s r)\ and p — 
<? a po exp[GM I \(? s r)], where c s is the isothermal sound speed, 
po is the matter density at the infinity (approximately equal to 
that on the external boundary). To ensure the gradual start- 
up (smooth initial conditions) we start to rotate star with very 
small angular velocity corresponding to the corotation radius 
at the external boundary R out , r cor = R out (GM /fi^) 1/3 and 
later gradually increase the angular velocity form that up to 
much larger values during 10 periods of rotation of the disk 
at r = 1. The (final) angular velocity of the star was differ- 
ent in different runs ranging from very small (corresponding 
to r cor = 10) up to very large (corresponding to r cor = 1) 
with most of runs done at intermediate value r cor — 3 (an- 
gular velocity of the star in dimensionless units fi» w 0.19). 
Spinning-up of the star leads to spinning-up of the magneto- 
sphere and the low-density corona. We checked that rotating 
dipole brings the corona to corotation in a very short time 
(with a speed corresponding to propagation of a torsional 
Alfven waves through the corona) . 

We divided the external boundary to the region of 
corona, < Ms k and the region of the disk, > 0disk- In the 
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Figure 2. Figure shows formation and evolution of conical winds with time. The background shows logarithm of the absolute value of the poloidal 
matter flux, lines - are selected magnetic field lines. 
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Figure 3. Different values shown for one moment of time, T = 500 in the whole region (top panels) and near the star (bottom panels). The 
background shows (from left to right): the logarithm of the absolute value of the poloidal matter flux, pv p , angular velocity, f2, poloidal current, 
J p = rB^, and entropy S. Lines are selected magnetic field lines. Thick line is a magnetic field line which is used for analysis of flow in the 
direction of conical winds. Vectors are velocity vectors. An arrow shows that the distance s along the line used in the subsequent plots starts at 
the disk. 
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Figure 4. Left panel: Forces along the part of the magnetic field line shown on the Fig. [3] starting from the disk up to the point when the 
line curves towards the star. Labels are: G-gravitational, C -centrifugal, M-magnetic, P-pressure gradient, and R - resulting forces. Right panel: 
velocities projected to the same line: vtot - total, vp - poloidal, - azimuthal, c s - sound speed, va - Alfven, v 3m - slow magnetosonic, Vf m - 
fast magnetosonic, v esc - escape velocities. 



corona we take "free" boundary conditions for all variables, 
but do not permit matter to flow back to the simulation re- 
gion from the boundary. In the disk region, we fix the density 
at some value, p = p d , and establish slightly sub-Keplerian 
velocity, Q,d = kQ (fd), where k — y/1 — 0.003 so that matter 
flows into the simulation region through the boundary. When 
the matter enters the simulation region it continues to flow 
inward due to the viscosity incorporated in the code. We use 
an a— type viscosity with viscosity coefficient proportional to 
a vis . Diffusivity is incorporated into the code with diffusiv- 
ity coefficient proportional to coefficient ctdif. Inflowing mat- 
ter has zero magnetic flux. The boundary conditions on the 
equatorial plane and on the rotation axis are symmetric and 
antisymmetric (see e.g. Ustyugova et al. 2006). 



of the star, R — 2R t and the surface magnetic field B t is 
different for different types of stars. The reference velocity 
is vq = (GMq/ Rq) 1 / 2 , the time-scale to = Rq/vq, and the 
angular velocity fi = l/t . We measure time in units of the 
rotational period with Keplerian speed at r = Ro : Po = 2wto 
and in the plots show dimensionless time T = t/Po. Other 
reference values: the reference density p and number density 
of particles no, the reference accretion rate Mo = povoR 2 , and 
the reference angular momentum flux Lo = MoVoRo- Table 1 
shows examples of reference variables for different stars. 

We solved the MHD equations for normalized variables: 
p = p/po, v = v/vo, B = B/Bq, etc. The plots show the 
normalized variables (with tilda's implicit). Examples in di- 
mensional variables are given for young stars. 



2.2. Reference Units. The MHD equations were solved in 
dimensionless form so that results can be applicable to stars 
of different scales, such as to classical T Tauri stars (CTTSs), 
to accreting white dwarfs, and to neutron stars. We chose a 
reference mass Mo = M» to be typical mass for different types 
of stars (see Table) . The reference radius is twice the radius 



3 PHYSICS OF CONICAL WINDS 

Multiple runs were done in the axisymmetric case for a wide 
variety of parameters. We chose one set of parameters as our 
main case and describe it in detail. We observed that the con- 
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Figure 5. Left panels: matter fluxes to the star M a tar, to conical winds M m i nd and through the disk Mdisk- Right panels: angular momentum 
fluxes carried by the disk, Ldi s k> to the star carried by matter, L ms tar, and by magnetic field, Lj star , and to conical winds, L mmin ds and 
Lf winds- Bottom panels show fluxes during a part of the simulation time. 
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Figure 6. Left panel: the background, vectors and streamlines show angular momentum carried by magnetic field. Middle panel: background 
and arrows show angular momentum carried by matter, streamlines show angular momentum carried by the field. Right panel: the background 
shows angular momentum caried by viscosity. 



dition ct vis > a d %s is necessary for conical winds to form. 
That is, the magnetic Prandtl number of the turbulence, Pr — 
OL V is/a.dif > 1. In addition, the conical winds are stronger if 
the diffusivity is not very small, that is if adif > 0.03. Below 
we show results for a vis — 0.3 and adif = 0.1 which we term 
"the main case". 



For the main case we chose densities in the disk and 
corona: pd = 10, p~ c = 0.001, corotation radius f cor — 3. 
Results can be applied to different types of stars. However ,as 
an example, we show results in dimensional units for typi- 
cal CTTS, where the reference values are taken from Table 1. 
For example, for a CTTS with a P» = 5.4 days, the unit of 
time used in figures below is Po = 1-04 days (see Table, CTTS 
column) . 



3.1 Matter flow, velocities and forces 

The simulations show that when the disk comes close to the 
star and bunches the poloidal field lines into an X-type con- 
figuration then conical outflows start and continue for a long 
time, as long as the simulations were run. The Fig. [l] shows a 
typical view of the matter flow where the matter flux is shown 
as a background. Matter flows to the cone into a relatively 
narrow shell and the cone has a half-opening angle, 9 — 30° — 
40°. The process is non-stationary due to episodes of inflation 
and reconnection of the magnetic field lines (see animation at 
http://www.astro.cornell.edu/~romanova/conical.htm One 
can see that the disk comes quite close to the star. Test simu- 
lations (which require much longer runs) have been done for 
an inner boundary smaller by a factor of two. This showed 
that the conical winds are similar, while the magnetosphere 
is larger. The poloidal velocity vectors show that matter is ac- 
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Figure 7. Density and velocity distributions in conical winds for CTTS. Here we converted dimensionless values obtained from simulations to 
dimensional values for typical T Tauri star taken from the Table. 
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Figure 8. Two components of winds from slowly rotating star are labeled. 



0.15 



AU 



celerated and forms a hollow conical wind. At the same time 
low-density matter is accelerated up to high velocities (hun- 
dreds of km/s in CTTSs case) along the inner surface of the 
cone (i.e., closer to the axis). This high-velocity component 
may be important in explanation of some highly blue-shifted 
spectral lines which form near CTTSs (Edward et al. 2003; 
2006). Matter which is accelerated in this region may come 
from the star and/or be captured from the main accretion 
flow. 

The Fig. [2] shows that the conical winds form and later 
become self-supported for a long time. Figure shows that all 
magnetic flux which threaded the disk at T = is bunched 
later by the disk into X-type configuration. The inward mo- 
tion of the disk is slow: it took about T — 100 rotational 



periods Po (at r = 1) to reach the vicinity of the star. The 
outflows started at time T £s 150 and continued for many 
hundreds of rotations. The magnetic field in the conical winds 
reconnects frequently and matter is expelled to conical winds 
in blobs with time-scale about (5 — 6)Po (see animation at 
http://www.astro.cornell.edu/~romanova/conical.htm ) . 

The Fig. [I] shows conical outflows at one moment of time 
T = 500. The figure also shows velocity vectors of the flow. 
Most of the matter is accelerated along conical winds. Some 
low-density matter is accelerated to much higher velocities on 
the inner side of the cone. 

The Fig. [3] shows the distribution of different parameters 
as a background at time T = 500 in the whole region (top) 
and in the inner part of the region (bottom) . The middle pan- 
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els show the distribution of angular velocity Q and lines of 
magnetic flux. One can see that the inner region of the closed 
magnetosphere rotates with the angular velocity of the star 
(1.2 < r < 2), while the outer region of the corona above the 
disk (at r > 2) rotates with the angular velocity of the disk. 
The field lines which start at the disk go through the regions of 
decreasing angular velocity and thus they are strongly wound 
up owing to the difference in the angular rotation rates. This 
leads to strong poloidal current rB,p above the disk (see next 
panels to the right) and the magnetic force associated with 
gradient of magnetic pressure, F m ~ —grad{(rB<f,) 2 ] appears. 
This is the main force driving matter to conical winds (as pro- 
posed by Lovelace et al. 1991). The right-hand panels (with 
entropy as a background) show that matter flowing to conical 
winds is cold. 

The Fig. [4] shows the projection of the different forces 
onto a field line (shown in bold on the bottom panels of Fig. 
[3]). This is one of the field lines which crosses the region where 
matter is accelerated to conical winds so that we calculate 
projection of forces to this line starting from the disk and ap- 
proximately half-way to the star (before it bends towards the 
star). The Fig. [4] shows that both, centrifugal C and gravita- 
tional G forces are large and approximately compensate each 
other. The centrifugal force is smaller, because a star rotates 
slowly, so that the sum G + C is negative. The pressure gradi- 
ent force P is positive inside and near the disk, but it is only 
strong enough to compensate negative value of G + C. The 
magnetic force M is the main one acting in the positive direc- 
tion and accelerating matter to conical winds. The gradient of 
this force is directed approximately in the z— direction, and 
matter is pushed vertically up. However, it is stopped by the 
magnetic flux of the dipole and is redirected into a cone. This 
is why the hollow conical flow is relatively thin and has a rel- 
atively small half-opening angle. If the centrifugal force were 
dominate (as in X- winds, Shu et al. 1994) , then the cone 
would have larger opening angle, and it would not be so nar- 
row, because the centrifugal force has horizontal direction. In 
addition, the vertical direction of the magnetic force leads to 
frequent forced reconnection of the magnetic flux threading 
the conical flows. 

The right-hand panels of Fig. [4] show the velocity vari- 
ation along the same same magnetic field line. One can see 
that the total velocity v to t is almost solely determined by the 
azimuthal rotation of the flow in the conical wind. Matter 
is accelerated in the poloidal direction from ~ c s (near the 
disk) up to slow magnetosonic v sm , Alfven va and fast mag- 
netosonic vf m velocities. The azimuthal component of veloc- 
ity v$ is much larger than poloidal component, > 3v p . 
The projection of the poloidal velocity to the chosen field line 
decreases at larger distances s because the field line curves 
towards the star. In reality, matter is accelerated more, up to 
~ v w 0.5 in dimensionless units. Note that the sound speed 
c a is small so that the pressure gradient is not significant in 
driving the conical winds. 

3.2 Matter and angular momentum fluxes 

We calculated the matter fluxes to the star's surface (inte- 
grated over R = 1), to the conical wind (calculated at R = 6), 
and the matter flux through the disk at R = 6. The Fig. [5] (left 
panels) show matter fluxes. One can see that most of matter 
incoming in the disk Mdisk accretes to the star, M sta r, while a 



smaller part (1/3 in the main case) flows to the conical wind. 
The bottom panel shows matter flux curves for a small inter- 
val of time for better temporal resolution. We also calculated 
the angular momentum fluxes to the star, to the conical wind, 
and through the disk. 

We calculated the flux of angular momentum through the 
surface of the star and through the surface R = 6. In both 
cases the flux was calculated as 



dS ■ r sin 9 



Ait 



v t pr sin^Vw 



where dS is the surface element directed outward to the re- 
gion. The first term on the right-hand-side L gives the trans- 
port of angular momentum by the matter; the second term is 
magnetic field contribution; and the third term is the trans- 
port due to viscous stress. We also calculated the density of 
the angular momentum flux and show it as a background in 
the Fig. [6] 

The Fig. [5] (right panels) shows that the largest flux is 
the inward flux carried by the disk. Only small part (^ 1/3) 
of this flux is carried away by conical winds. Another small 
part («i 1/4) is transported to the (slowly rotating) star 
mainly through magnetic interaction with much smaller por- 
tion transported by the matter. The rest of the flux is trans- 
ported back by viscous stresses. 

The Fig. [6] shows angular momentum flow. The left panel 
shows that angular momentum to the star is transported by 
magnetic field. Most of transport occurs through magnetic 
field lines in the area of the funnel stream, where most of 
matter flows (see also Romanova et al. 2002; Bessolaz et al. 
2007). One can see that the magnetic stress is also high at the 
base of the conical outflows which means that it is responsi- 
ble for transport of the angular momentum at the base of the 
conical winds. Middle panel shows that later, above the disk, 
angular momentum is carried by matter. Right panel shows 
angular momentum carried by viscous stress. 



3.3 Application to CTTSs 

CTTSs are strongly variable on different time-scales including 
a multi-year scale (Herbst et al. 2004; Grankin et al. 2007). 
This is connected with variation of the accretion rate through 
the disk which may lead to the enhancement of outflows (e.g., 
Cabrit et al. 1990). Simulations have shown that the bunch- 
ing of field lines by the new matter after period of the low 
accretion may lead to quite long outburst of matter to coni- 
cal winds and may be the reason for formation of micro-jets 
in the CTTSs. If CTTS is in a binary system, then an accretion 
rate may be episodically enhanced due to interaction with the 
secondary star. Events of fast, implosive accretion are possi- 
ble due to thermal instability or global magnetic instability, 
where the accretion rate is enhanced due to the formation of 
disk winds (Lovelace et al. 1994). Fig. [7] shows an example 
of conical winds in case of CTTSs. Some observations confirm 
a hollow cone shape of the outflows from CTTS. For exam- 
ple comparison of possible geometries of outflows observed 
in Hp line in RW Aurigae led to the conclusion that a conical 
shaped wind with half-opening angle 30 — 40° and a narrow 
annulus gives the best match to the observations of this line 
(from Alencar et al. 2005; see Fig. [9]!. 

Two components of the outflow are observed (see Fig. [8]) 
- a slow component, associated with the conical wind, and a 
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fast component associated with accelerated low-density mat- 
ter closer to the rotating axis. This high- velocity component 
may be important in understanding some highly blue-shifted 
spectral lines which form near CTTSs. Matter which is accel- 
erated in this region may come from the star or it may be cap- 
tured from the main accretion flow. Observations show that 
a significant number of CTTS show signs of outflows in spec- 
tral lines, in particular in He I (Edwards et al. 2006; Kwan, 
Edwards, & Fischer 2007). 

If CTTS typically in rotational equilibrium state (e.g. 
Long et al. 2005), then the stars will oscillate between periods 
of accretion and periods of propeller-driven outflows. In the 
propeller stage the high velocity component of the outflow is 
even stronger (see §5). 



3.4 Periods of enhanced accretion and outflows in 
EXors 

A period of enhanced accretion may lead to outbursts in 
EXors, where the accretion rate increases up to 10~ 6 — 
10 _5 M Q //yr and strong outflows are observed. Brittain et al. 
(2007) reported on the outflow of warm gas from the inner 
disk around EXor VI 647 observed in the blue absorption of 
the CO line during the decline of the EXor activity. He con- 
cluded that this outflow is a continuation of activity associ- 
ated with early enhanced accretion and bunching of magnetic 
field lines (see Fig. |10p . 

Our simulations are directly applicable to such a situa- 
tion: new matter comes after period of low accretion rate in 
CTTS and magnetic field lines are bunched close to the star. 
The accretion rates are higher than in CTTSs so that all di- 
mensional parameters will correspond to stronger outflows. 
In our main example of a CTTS, the disk stops at R m — 2AR*. 
In EXors, we take the radius of a star at the Fig. [7] equal to the 
inner boundary, so that the disk stops much closer to the star, 
R m = 1.27?* . Then all velocities are a factor 1.4 higher and 
densities a factor of 32 higher (compared to Fig. [7]), and the 
matter flux in Fig. [5] is a factor of 1 1 higher than in the main 
example relevant to CTTSs. 



3.5 Comparison of the model with outflows from the 
Circinus X-l - the neutron star hosting binary 

Circinus X-l represents a case where jets are seen from the 
vicinity of the accreting neutron star. The system is unusual 
in the sense that Type I X-ray bursts as well as twin-peak X- 
ray QPOs are observed. The neutron star is estimated to have 
a weak magnetic field (Boutloukos et al. 2006). The binary 
system in this case has a high eccentricity (e ~ 0.4 — 0.9) and 
thus has periods of low and high accretion rates (e.g., Mur- 
din 1980). Two-component outflows are observed. The radio 
observations show a non-stationary jet on both arcmin and 
arcsecond scales, which appears to have a small opening an- 
gle. Recent spectroscopic observations in optics (Jonker et al. 
2007) and in X-ray band (Iaria et al. 2008; Schulz et al. 2008) 
have shown that outflows have a conical structure with a half- 
opening angle, about 30°. Different explanations are possible 
for such a conical structure, such as precession of a jet (Iaria 
et al. 2008). However, this appears less likely because the axis 
of the jet did not change in 10 years (Tudose et al. 2008). In 
our understanding this neutron star is a good candidate for 
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Figure 9. Modeling of the Hg line in RW Aurigae led to the conclu- 
sion that a conical shaped wind with half-opening angle 30 — 40° and 
a narrow annulus gives the best match to the observations of this line 
(from Alencar et al. 2005). 



conical outflows, because (1) it has episodes of very low and 
very high accretion rates, which is very favorable for forma- 
tion of conical outflows, (2) a neutron star has only weak 
magnetic field which can be strongly compressed almost to 
the stellar surface by the disk which is favorable for conical 
outflows. This also means that velocity of outflows are of the 
order of the Keplerian velocity at the surface of the neutron 
star. 



3.6 Application to black hole hosting systems 

Jets and winds are observed from accreting black holes in- 
cluding both, stellar-mass black holes and black holes in 
galactic nuclei. The correlation between enhanced accretion 
rate and outflows has been extensively discussed and observa- 
tional data are in favor of this correlation. On the other hand 
the possibility of magnetic flux accumulation in the inner disk 
around the black hole and the enhancement it gives to jets 
has been discussed and observed in numerical simulations 
(Lovelace et al. 1994; Meier 2005; Igumenshchev 2008). In 
this situation periods of enhancement of the accretion rate 
may also lead to bunching of the field lines into an X-type 
configuration (near the inner radius of the disk) and to for- 
mation of conical winds, because the mechanism of conical 
winds does not require the dipole or any other special config- 
uration. Similarly, the magnetic field of a star may be more 
complex than a dipole one (see e.g., Donati et al. 2006; Mo- 
hanty and Shu 2008). 



4 OUTFLOWS IN THE "PROPELLER" REGIME 

In case of accreting magnetized stars outflows may be asso- 
ciated with "propeller" stage of evolution (e.g., Illarionov & 
Sunyaev 1975; Alpar & Shaham 1985; Lovelace et al. 1999) 
if a star rotates sufficiently fast. At propeller stage fast ro- 
tating magnetosphere transfers its angular momentum to the 
disk matter (e.g. Lovelace et al. 1999). Simulations show that 
most of matter flows from the inner parts of the disk to widely 
open cones (Romanova et al. 2005; Ustyugova et al. 2006). 
The process is quasi-periodic, or episodic (with no definite 
period) which is connected with opening and closing of mag- 
netic field lines connecting a star and the disk. Such quasi- 
periodic reconstruction of the magnetosphere due to inflation 
and reconnection has been predicted theoretically (Aly & Kui- 
jpers 1990, Uzdensky, Litwin & Konigl 2003) and has been 
observed in a number of axisymmetrtic simulations (Hirose 
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Figure 10. Schematic model of an Exor VI 647 Ori. During the outburst the accretion rate is enhanced so that the magnetospheric radius R m 
decreases and the magnetic field lines were bunched (A) . This results in a fast, hot outflow. As the accretion rate decreases, the disk moves outward 
and this results in a slower, cooler CO outflow (B). Further decrease in the accretion rate leads to a quiescence state where the production of 
warm outflows stops (C). From Brittain et al. (2007) 




Figure 11. Outflows in the propeller regime. The background shows matter flux, lines are selected field lines, arrows are proportional to velocity. 
Labels show total velocity and density at sample points. 



et al. 1997; Goodson et al. 1997, 1999; Matt et al. 2002; 
Romanova et al. 2002; von Rekowski & Brandenburg 2004). 
Multiple events of inflation and reconnection (hundreds of 
outbursts) were observed in simulations of the propeller stage 
(Romanova et al. 2005; Ustyugova et al. 2006; see animation 
at http://www.astro.cornell.edu/~romanova/propeller.htm ) . 

In the propeller regime the magnetosphere rotates faster 
than inner region of the disk. This occurs if the co-rotation 
radius R cr = (GM/fi*) 1 ' 3 is smaller than magnetospheric 
radius R m (e.g., Lovelace et al. 1999). Young stars are ex- 
pected to be in the propeller regime in two situations: (1) At 
the early stages of evolution (say, at T < 10 6 years), when 
the star formed but did not have time to spin-down, and (2) 
at later stages of evolution, such as at CTTS stage, when the 
star is expected to be on average in the rotational equilibrium 



state (e.g., Long et al. 2005) but variation of the accretion 
rate leads to variation of Rm around R cr , where R cr < R m 
is possible. We performed axisymmetric simulations of accre- 
tion to a star in the propeller regime, taking a star with the 
same parameters as in case of conical winds, but with period 
P* = 1 day (Romanova et al. 2005; Ustyugova et al. 2006). 
We chose a v — 0.3 and a d = 0.1 and thus bunched the field 
lines to the X-type configuration We observed that in addition 
to conical wind there is a fast axial jet (see Fig. \\\\ so that 
the outflow has two components (see Fig. |12| l . The conical 
wind in this case is much more powerful - it carries most of 
the disk matter away. The axial jet carries less mass, but it 
is accelerated to high velocities. Acceleration occurs due to 
the magnetic pressure of the "magnetic tower" which forms 
above the star as a result of winding of magnetic field lines 
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of the star. Outbursts to conical winds occur sporadically with 
a long time-scale interval (see Fig. |13[ l which is connected 
with the long time-scale interval of accumulation and diffu- 
sion of the disk matter through the magnetosphere of the star 
(see also Goodson et al. 1997; Fendt 2008). These propeller 
outflows were obtained in conditions favorable for such a pro- 
cess: when the star rotated fast and an X-type configuration 
developed. Future simulations should be done for the case of 
propeller-driven outflows from slower rotating CTTS. Colli- 
mation of conical winds may occur at larger distances from 
the star for example, by disk winds (e.g., Konigl & Pudritz 
2000; Ferreira et al. 2006; Matsakos et al. 2008). 



5 3D SIMULATIONS OF CONICAL WINDS 

We performed exploratory simulations of conical winds in 
global 3D simulations. One of main questions is what is an 
expected direction of conical winds when magnetosphere of 
the star is misaligned relative to the rotational axis. We cho- 
sen a case when magnetosphere of a star is misaligned rel- 
ative to the rotational axis at an angle O = 30°. We took 
a cubed sphere code developed earlier (Koldoba et al. 2002, 
Romanova et al. 2003) and used a grid N r xN 2 = 120 x 51 2 in 
each of 6 blocks of the sphere. At such fine grid we were able 



to handle the case with the low-density corona p c = 0.001 
and took the density in the disk pd — 2 which is 5 times 
lower than in the axisymmetric case. However we also took 
smaller magnetic moment of the star, fx = 2 (compared to 
fb = 10 in the axisymmetric case), again, to save comput- 
ing time. We did not push a disk from too large distances (to 
save computing time) but from r — 5 which is far enough 
to insure an initial bunching of the magnetic field lines. Sub- 
sequently, the bunching had been supported by sufficiently 
high viscosity, a vis = 0.3. We do not have diffusivity in the 
3D code however at such a grid estimated numerical dif- 
fusivity at the disk-magnetosphere boundary is at the level 
a d ~ 0.01 — 0.02 so that conditions for conical winds are sat- 
isfied, that is Pr m > 1. 

Simulations have shown that accreting matter bunched 
field lines and that some matter went to conical winds. Fig. 
^] shows that conical winds are geometrically symmetric rel- 
ative to the rotational axis of the star Q, in the sense that an 
angle between different parts of the conical winds and fl axis 
is approximately the same. However the density distribution 
in the wind shows a spiral structures in the form of density 
amplification, which rotates with an angular velocity of the 
star, fit,, and represents a one-arm spiral from each side of 
the outflow. 

Note that at a vis — 0.3 the disk-magnetosphere interac- 
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Figure 14. Projections show formation of conical winds obtained in 3D MHD simulations. Left panels show an XZ slice in the distribution of 
matter flux (background) and magnetic field lines in different moments of time. Right panels show corresponding YZ slices. Arrows show the 
direction of the magnetic moment of the star /j and angular velocity of rotation Q. 



tion has tendency to be unstable relative to the interchange 
instability (Romanova, Kulkarni & Lovelace 2008; Kulkarni & 
Romanova 2008) and we do observe partial accretion through 
instabilities in addition to funnel stream accretion which dom- 
inates at pretty high inclination of the dipole = 30° (Kulka- 
rni & Romanova 2009). One can see that conical winds start 
blowing from those parts of the disk which are quite far 
away from the disk-magnetosphere boundary. We can guess 
that again, the accumulation of the component above the 
disk drives matter outward to the conical winds. However 
we should be caution in giving global conclusion about non- 
dependence of conical winds from instabilities. Penetration of 
the inner disk matter through the magnetosphere through in- 
stabilities decreases the density at the inner edge of the disk 
and ability of the disk matter to strongly compress the mag- 
netosphere which is a favorable condition for strong conical 
winds. This possible interrelation between instabilities and 
conical outflows should be checked in the future 3D simu- 
lations in greater detail. 



6 INVESTIGATION OF PARAMETER RANGE 

To investigate the dependence on different parameters we 
took the main case and varied one parameter at a time. The 
main effort has been directed at understanding the depen- 
dence on the effective viscosity and the effective magnetic 



diffusivity. We performed two sets of runs: (1) at fixed dif- 
fusivity and different viscosity coefficients, and (2) at fixed 
viscosity and different diffusivity coefficients. 

Dependence on viscosity at fixed diffusivity. We fixed the 
diffusivity at a di f — 0.1 and varied the viscosity coefficient 
in the range a viB = 0.01 — 1. We observed that for small vis- 
cosity, a vis < 0.1, the magnetic field of the dipole diffuses 
through the inner regions of the incoming disk and an X-type 
configuration does not form. No conical outflows appear in 
this case. We conclude that formation of conical winds re- 
quires a viB > adif, that is, Pr m > 1. Next we increased 
a vis and observed that the X-type configuration formed and 
conical outflows were generated. We observed that the ac- 
cretion rate to the star strongly increases with a viB , while 
outflow rate to the conical winds increases but only slowly. 
For a V i S — 0.1, the matter fluxes to the star and to the wind 
are small and approximately equal. For a vis — 0.3&0.4, the 
wind carries about 30% and 20% of mass correspondingly). 
Angular momentum carried to the star also strongly increases 
with a vis . In all cases the star spins-up, because the magne- 
tospheric radius, r m w 1.2 is smaller than corotation radius 
which is r CO T = 3. That is, incoming matter brings positive an- 
gular momentum to the star. Conical outflows carry angular 
momentum away from the disk. We should remember how- 
ever that it is only a small part of the total angular momentum 
of the disk as shown above. For a vis — 0.3 and 0.4, matter 
outflows to the conical winds in oscillatory fashion. We often 
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observe that quasi-periodic oscillations of the accretion occur 

for Pr m w 3 - 4. 

Dependence on diffusivity at fixed viscosity. In the next set 
of runs we fixed the viscosity at a vi3 = 0.3 and varied the dif- 
fusivity: from O.Oliol. Again, no conical winds were formed 
for ctdif > a vis . At relatively high diffusivity, adif = 0.1, 0.3, 
about 30% of the incoming matter flows to the conical winds. 
For ctdif = 0.1 (Pr m = 3) conical outflows oscillate, while 
at adif = 0.3 (Pr m = 1) no oscillations are observed. For 
very small diffusivity, adif = 0.01, the conical outflows also 
form with slightly smaller matter flux to the winds. However, 
the accretion rate to the star is larger. This is somewhat un- 
expected because one would expect that at small a d i / matter 
of the disk would be trapped by strong magnetic flux accu- 
mulated at the inner edge of the disk. The possible explana- 
tion may be that penetration of matter of the disk through 
the disk-magnetosphere boundary is determined by the lo- 
cal Reynolds number Re m = A.rv r /r) m , which is determined 
not only by the diffusivity coefficient j) m ~ adif, but also by 
the velocity of the flow, v r which become very small at the 
disk-magnetosphere boundary, and also by the size A r of the 
boundary, which also become very small at high Prandtl num- 
bers. Angular momentum fluxes to the star and to the winds 
are approximately the same, excluding the case adif = 0.01 
where the flux to the star is larger. This approximate equality 
does not mean that the conical winds carry angular momen- 
tum flux out of the star. Instead, they carry angular momen- 
tum from the inner regions of the disk. The remaining angular 
momentum of the disk goes to the spinning-up the star and 
part of the angular momentum of the disk flows outward in 
the disk due to the viscous stress. 

Variation of other parameters. We varied the period of the 
star (corotation radius) taking r cor = 5, 10 for slowly ro- 
tating stars and r cor = 1.5, 2 - for faster rotating stars. In 
the case of slowly rotating stars conical winds form and the 
outflow rate to the winds is similar to one in the main case 
for r cor — 3, though accretion rate to the star is somewhat 
larger. For faster rotating star, r cor = 2, the amplitude of vari- 
ability increased and matter flux to outflows increased up to 
50%— 70% compared to the accretion rate to the star. For even 
faster rotation, r cor — 1.5, the accretion rate decreased by a 
factor of 5 compared to the main case, while the outflow to 
conical winds increased by a factor of 2. Thus most of matter 
is ejected to conical winds. This situation become close to the 
propeller regime, where the corotation radius becomes close 
to the magnetospheric radius, r m « 1.2 where most of the 
incoming matter may be ejected to the outflows (Lovelace et 
al. 1999; Romanova et al. 2005; Ustyugova et al. 2006). 

Outflow of matter to winds occurs if the corona is not 
very dense so that outflowing matter of the winds does not 
lose its energy while propagating through the corona. In the 
main simulation runs an initial density of the corona is 10~ 4 
times the disk density (p d = 10 versus p c — 10~ 3 ). To test 
the dependence on the coronal density we decreased its den- 
sity by a factor of 3 and first chose smaller (compared to the 
main case) transport coefficients, a vis = 0.1, a di f = 0.03, 
hoping to enhance outflows to the winds. Simulations have 
shown that matter fluxes to the star and to the winds are not 
appreciably different from the main case. We conclude that 
the coronal density used in the main case is sufficiently small 
as to not suppress the outflows. 



Summarizing the above subsections we conclude that: 

(1) Strong conical outflows appear at a wide range of pa- 
rameters if both the magnetic viscosity and diffusivity are not 
very small, a v > 0.1 and ay > 0.01. Outflows are most pow- 
erful when the viscosity is a few times larger than diffusivity. 
For these parameters the viscosity is high enough to drive disk 
matter inward to the region of stronger magnetic field, while 
the diffusivity is high enough to ensure the penetration of disk 
matter through the bunched field lines. 

(2) Weak outflows are observed when the viscosity is large 
and the diffusivity is very small. In this case the field lines are 
bunched into the X-type configuration. However, the diffusiv- 
ity is not high enough to ensure loading of the disk matter to 
the field lines of the magnetosphere. 

(3) Slow accretion, no outflows occurs when the viscosity 
is small, a v < 0.1 for any value of the diffusivity a d . In this 
case the disk is stopped at larger distances from the star by 
the strong magnetic field of the star. Only very week outflows 
or no outflows are observed in this case. 



7 CONCLUSIONS 

We discovered a new type of outflows - hollow conical winds 
- in MHD simulations of disk accretion to a rotating magne- 
tized star. The conical winds occur under conditions where 
the poloidal magnetic field is bunched into an X-type config- 
uration. 

In some respects these winds are similar to the X-winds 
proposed by Shu and collaborators (e.g., Shu et al. 1994): 

1. They both require bunching of the field lines; 

2. They both have high rotation of the order of Keplerian rota- 
tion at the base of outflow, and gradual poloidal acceleration; 

3. They both are driven by the magnetic force (see also 
Lovelace et al. 1991). 

However, there are a number of important differences: 

1. Conical winds flow in a thin shell (they are hollow), while 
X-winds flow at different angles below the "dead zone"; 

2. Conical winds form around stars of any rotation rate in- 
cluding slow rotation, and do not require the fine tuning of 
angular velocity of the inner disk to that of magnetosphere; 

3. Conical winds are non-stationary: the magnetic field con- 
stantly inflates and reconnects; 

4. Conical winds carry away part of the angular momentum 
of the inner disk and are not responsible for spinning-down 
the star, while X-winds are predicted to take away angular 
momentum from the star; 

5. In the conical winds there is a fast component of the flow 
along field lines threading the star. 

6. In the propeller regime (where the corotation radius is 
comparable or smaller than magnetospheric radius) a new 
component appears: a strong magnetically dominated jet 
along the open field lines of the star. In this component low 
density matter is accelerated up to super-Keplerian veloci- 
ties. The star loses its angular momentum mainly due to the 
twisted magnetic field of this component. 

Some of these differences, such as non-stationarity of conical 
winds is connected with natural restrictions of the stationary 
model of X-winds. 
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Conical winds can explain the hollow conical shape of 
outflows near young stars of different type (CTTSs, EXors, 
Type I objects) which have been recently resolved. The coni- 
cal winds may arise from any type of the disk-accreting mag- 
netized stars, including accreting young stars, white dwarfs, 
neutron stars and possibly black holes as mentioned. An ex- 
act configuration of the magnetic flux is not important: only 
the bunching/inclination and inflation of the flux are impor- 
tant. This makes this mechanism universal and it may work 
in different situations. It is often observed that enhancement 
of the accretion rates leads to enhanced outflows. Mechanism 
of conical winds gives natural explanation of such events ob- 
served in a wide variety of astrophysical objects. 

It is still not known what are the values of transport coef- 
ficients in real accretion disks. MRI simulations have shown 
that magneto-rotational mechanism of turbulent transport 
may give values of a in the range: a vis — 1CP 2 — 0.4 (e.g., 
Stone et al. 2000) and corresponding diffusivity values. If in 
real disks transport coefficients are large, a ~ 0.1 then strong 
outflows are expected during periods when a vis > otdi / with 
steady flow to the wind or with short-period variability like 
in Fig. [5] If in the opposite case the transport coefficients 
are always small, say a viB = otdif = 0.01, then conical out- 
flows are expected only during periods of enhanced accre- 
tion when accretion rate is enhanced due to e.g. some insta- 
bilities or other mechanism of angular momentum transport, 
e.g., Rossby waves (Lovelace et al. 1999; Li et al. 2000). Then 
outflows will be characterized by high-amplitude outbursts 
to conical winds with a large time-interval between outbursts 
determined by the small diffusivity at the disk-magnetosphere 
boundary (similar to that in the propeller regime (e.g., Fig. 

LH. 

If a star is in the propeller regime, then clear two- 
component outflows are observed: 

(1) Conical-type winds which may carry away significant part 
of the disk mass and angular momentum, with velocities ap- 
proximately equal to Keplerian velocity at the base of the 
winds. 

(2) Magnetically dominated axial jet where low-density mat- 
ter is accelerated near the star to much higher velocities. The 
twisted magnetic field of this jet carries away most of the an- 
gular momentum of the star. 

The propeller regime may be responsible for the fast loss 
of angular momentum by young stars before their CTTS stage. 
It may also appear in stars which are in rotational equilibrium 
during periods of the accretion rate is lower. 

The propeller-driven axial jet has good collimation due to 
interaction with surrounding conical wind. The overall coni- 
cal wind may be collimated at larger distances by disk winds 
emanating from the disk at large radii (Blandford & Payne 
1982) or the pressure of the external medium (Lovelace et al. 
1991). Typical terminal velocity of outflows from the disk are 
of the order of the Keplerian velocity at the radius from which 
matter flows to the wind and is expected to be smaller than 
velocity in conical winds, though the total angular momen- 
tum maybe larger (e.g., Ferreira et al. 2006). 
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